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Abstract 


A kinetic model for the anode of the direct methanol fuel cell (DMFC) is presented. The model is based on the generally accepted dual site 
mechanism of methanol oxidation, in aqueous solution, on well characterized Pt—Ru catalyst and it can predict the performance of the electrode 
as a function of cell temperature, anode potential and methanol concentration. In addition the model also generates data regarding the surface 
coverage of significant adsorbates involved in methanol oxidation on the dual site catalyst. 

The analysis of the initial complex model confirms that a simplification in anode modelling can be made and some of the kinetic parameter can 
be reliably neglected. Based on this approach a fast and simplified three parameter model is derived from the same complex kinetic mechanism. The 
kinetic parameters of both models are estimated from experimental anode polarisation data from a 9 cm? DMFC operating with various methanol 
feed concentrations and temperatures. The models were developed in Lab VIEW and this has greatly simplified the simulation process, giving a 
model with ca. 85-95% fit on the experimental data. Depending on the computational speed available, and the desired complexity of problem at 


hand, either of the models can be used to give accurate model simulations for methanol fuel cell polarisations. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Low temperature liquid feed direct methanol fuel cells 
(DMEFCs) are promising power sources for portable, stationary 
and vehicular applications given the relatively compact system 
design and higher energy densities when compared to existing 
technologies [1-3]. However, obstacles, such as the relatively 
poor kinetics of methanol oxidation at the anode, crossover 
of methanol through the membrane and the subsequent mixed 
potential at the cathode still hinder their widespread commer- 
cialisation [1-5]. 

The kinetics of the methanol oxidation reaction are deter- 
mined by a complex mechanism involving adsorption of 
methanol on catalyst site followed by parallel and/or series reac- 
tions for electrochemical oxidation of methanol to CO2. Based 
on the pioneering work of Bagotzky et al. [6] several mechanisms 
have been proposed in the literature to describe the methanol 
oxidation process [3,7—13]. Modern in situ spectroscopy meth- 
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ods like FTIRS and mass spectrometry have shed more light 
on the surface intermediates formed during the reaction and 
have been employed to elaborate upon the mechanism under 
different operating conditions. The Pt—Ru alloy is considered 
to be the present state-of-the-art candidate for practical anode 
catalysis in the DMFC and the common consensus is that the 
electro-oxidation of methanol on Pt—Ru occurs via a dual site 
mechanism [1—3,9-11,14—20,22]. 

To understand the complex kinetics in this work a semi 
empirical model based up on a dual site kinetic mechanism is 
presented. The model highlights the limiting reaction based on 
the kinetic parameters and elucidates the surface coverage of 
intermediate species formed in the reaction. The kinetic param- 
eters deduced from the initial intricate mechanism are presented 
and further interpretation of these kinetic parameters has led to 
the formation of a simple kinetic expression which can be solved 
rapidly and used in real time DMFC simulations. 


2. Experimental 
In situ anode polarisation data was collected using an a 


single cell graphite (Ralph Coiden) fuel cell assembly etched 
with seven parallel channels of length 30mm x 1 mm x 2mm 
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Nomenclature 


a activity 

a/o atomic ratio 

C concentration (mol m~?) 

E potential (V) 

E, activation energy (kJ mol~!) 

F Faraday’s constant (96487 A s mol`!) 
current density (A m~?) 


J 

ki reaction rate constant for kinetic Eq. (1) 
k; reversible rate constant for kinetic Eq. (1) 
n number of electrons transferred 


r rate of reaction 
R gas constant (8.314 J mol`! K7!) 
T temperature (K) 


Greek letters 


r site density (mol m~?) 

Qj electrochemical transfer coefficient for kinetic 
Eq. (1) 

Bi transfer coefficient for kinetic Eq. (1) 

Oi surface coverage for intermediate i 

Subscripts 


ads adsorbed 

CO, Pt carbon monoxide on platinum sites 
OH, Pt hydroxyl ion on platinum site 

OH, Ru hydroxyl ion on ruthenium site 
H20 water 

M methanol 


creating an active surface area of 9cm?. An in-house manu- 
factured membrane electrode assembly (MEA) was employed 
to generate the necessary data [15,16]. The MEA was fab- 
ricated using 1mgp,;cm7? at the anode and cathode with 
Pt-Ru/C 60wt.% a/o 1:1 and Pt/C 60wt.% from Etek, respec- 
tively. Catalyst layers contained 20wt.% Nafion® ionomer and 
catalyst layers were applied upon a micro-porous layer of 
0.5mgKBcm~? (Ketjen Black) with a 20wt.% PTFE binder 
content deposited upon a TGP 090 20% WP (Toray) gas dif- 
fusion layer. The Electrodes were hot pressed to a Nafion® 117 
membrane at 140°C under a load of 50kgcm™? for 10 min 
and allowed to cool under pressure. The MEA was mounted 
inside the test cell and compressed to a torque of ca. 2 Nm. 
Conditioning over 24 h with intermittent polarisation resulted in 
reproducible performance. 

Anode data was collected by using the cathode as a pseudo- 
reference electrode, by feeding a 5%.vol. hydrogen in nitrogen 
stream and polarising the cell between OCP and 700 mV [16] 
using liquid methanol as fuel at the anode side. Data were col- 
lected potentiostatically and at a potentiodynamic sweep rate 
of 2mVs~“, with minimal difference observed. Conditions of 
temperature and concentration were varied between 30-90 °C 
and 0.254 mol dm~?, respectively, and the system were given 
ca. 2h to equilibrate before testing. 


3. Model development 
3.1. Reaction mechanism 


Itis widely accepted that during methanol oxidation, the most 
significant reactions are the adsorption of methanol on the cata- 
lyst site and the oxidation of CO. Hence, in this work, a simplified 
and general reaction mechanism derived from a more complex 
mechanism was selected [8,9,18]. The simplified mechanism 
using dual site approach comprises the following steps: 


3.1.1. On Pt sites 


CH30H > CH30Haas (1) 
1 
CH3 OH ads <=> COaas, pi + 4H" + 467 (2) 
2 
k3 Pi + = 
H20<— OHaas, Pt +H' +e (3a) 
3,Pt 


3.1.2. On Ru sites 
k3 Ru 
H2045 OHads,Ru + Ht + e7 (3b) 


3,Ru 


3.1.3. Surface reaction 
k. = 
COaas,Pt + (OHads, Pt + OH ads, Ru) > CO, +H* +e (4) 
‘4 


It is assumed that reaction (3b) principally occurs on ruthe- 
nium (Ru) and reaction (1) occurs principally on platinum (Pt), 
as it is known that Ru is a poor electrocatalyst for methanol 
oxidation. Generally it is thought that the rate controlling step 
is reaction step (4) between the adsorbed species COads,pt and 
OHaas, which in turn depends on reaction steps (1-3) for the 
formation of adsorbed species. Thus, on the basis of this the rate 
of reaction can be written as: 
r4 = k48ou,PrOco, pi el POPE RT) 


+ k400H, Ru9co, pr et PYFED/RT) (5) 


3.2. Assumptions 


The assumptions used in the model are 


i. The electrodes are at steady state. 

ii. Isothermal and isobaric operation and the concentration 
gradient between the reacting site and the bulk fuel is 
negligible. 

iii. Butler—Volmer kinetics is valid for the charge transfer steps. 

iv. Only liquid phase is considered and CO is assumed to be 
dissolved in the solution. 

v. Methanol oxidation described by the kinetic expression 
Eqs. (1)-(4) follows the dual site mechanism thus surface 
coverage are independent of the adsorption rate of other 
species. 
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3.3. Analytical solution 


The rates of changes of surface coverage of different species 
(0;) with respect to time are: 


d 


dé 
r a = iby e((@2FE)/RT) _ k59co.Pt e(((l—B2)FE)/RT) 


— kgoun,Pr9co,Pt € 


temperatures and concentrations adjusting seven rate constants, 
which presents a vast number of permutations. The model was 
therefore constructed using Lab VIEW software from National 
Instruments [21]. An alternative approach in which adsorption 


dôM = 
Pe = kıCm(1 — 60n,Pt — Oco,Pt — Om, Pt) — ky OM, pi + kb Aco, pr eLF2FEY RD) Oy pi eO2FE/ RT) (6) 


((1—f4)FE)/RT) ((—B4)FE)/RT) (7) 


— k400H,Ru bco, Pt € 


dôoH, Pt m Z = 
r 7. k3 pray,o(1 — 60n,Pt — 8co,Pt Ou pe PFE RE) k3 p,Oou,Pt e((CP3)FE)/RT) _ k49oun,Pt9co,Pt e((C— Ba FE)/RT) 
(8) 
don, R: _ _ _ 
r 7 a k3 Ruam o(1 _ bon Re! B3)FE)/RT) _ kh RuĝOH,Ru e B3)FE)/RT) _ k400H, Rub CO, Pt eC B4)FE)/RT) (9) 


Eq. (3a), the dissociation of water and adsorption of hydroxyl 
species upon platinum is known to be limited to relatively high 
potentials. Cyclic voltammetry and cell anode polarisation data 
have confirmed that OH-groups are preferentially formed on 
ruthenium at fairly low potentials (ca. 0.3 V versus dynamic 
hydrogen reference electrode (DHE)), whereas much higher 
potentials (ca. 0.9 V versus DHE) are required up on platinum 
[2,3,10,14,19,20]. On this basis, it can be assumed that adsorp- 
tion of hydroxyl] ions on Pt sites is negligible (0op,pt ~ 0) and 
can be neglected from all the above equations. Thus, Eq. (5) is 
simplified to: 


r4 = kaon Ru8co, pr el FYPP RD) (10) 

The case when adsorption of hydroxyl species also occurs on 
Pt, reaction (3a), is considered in the Appendix A along with 
the solution. This modelling procedure will contain additional 
three kinetic parameters hence it will be more time consuming 
and complex compared to the detailed model in manuscript. 

At steady state surface coverage does not vary with time 
and the derivative becomes zero. Thus, with ay,o = 1, and 
Aon Pt © 0, the steady state solution of Eqs. (6)—-(9) gives: 


kıCm — Om(ki Cm + kz e(@2FEV/RT) + ki) 11 
kı Cm — k, eC -AFE/RT) i 


0co,Pt = 


koOM e((@2FE)/RT) 


9co,Pt = ki, (OBERT) +. k40OH, Ru el BaFEY/ RT) (12) 


k3 pu 00 PIFEYRT) — k3 pyOou,Ru eC PERRI) 
a k3 RyOOH,Ru e(((—A3 )FE)/RT) 


kaon Ru eC BaFE)/RT) 


OCO, Pt = 
(13) 


The above simultaneous equations can be solved to give the 
surface coverage of 0OH,Ru and 6co,pt, and thus to find the overall 
rate of reaction (74). This can then be combined with Faradays 
law to give: 


z. j = n Fk480H, Rub co, pr FYFE RD) (14) 


The complexity in the solution of this problem arises as 
the quadratic must be solved at multiple potentials, at various 


of OH on Pt is also considered, is given in the Appendix A. It 
should be noted that addition of more reaction steps increases 
the number of parameters to be analysed and adds complexity 
in solving the non-linear equation [12,13,17]. The basis of semi 
empirical modelling is generally based on principal of parsi- 
mony and fast computational speed. A detailed model, as shown 
in the Appendix A, is complex and slow to compute. However, 
detailed models when supplement with transient electrochemi- 
cal measurements from modern tools can lead to the extraction 
of accurate kinetic parameters and limiting steps [3,20]. An addi- 
tional factor in the model used here is that the methanol oxidation 
is considered as a single four electron step. Introduction of single 
electron transfers involving more intermediates adds to the com- 
plexity and creates a large number of parameters which generally 
would not be amenable to estimation using the experimental data 
available. 


4. Results and discussion 
4.1. Model validation 


To date the trend in validation has been quite varied. Several 
authors verify their models using cell polarisation data alone 
without a reference electrode [1,2] whereas others have veri- 
fied their anode models against anode polarisation data using 
a hydrogen evolving cathode [4,5]. However, polarisation data 
alone, without any consideration of whether the model fits the 
Tafel slope for the data or whether the model can predict the acti- 
vation energy within a realistic range has rarely been reported. 
The model was deemed to provide a suitably accurate fit when 
the R? statistics for the actual and calculated curves was greater 
than 90%. These additional considerations, linked with the fact 
that coverage of the catalyst is restricted to values between 0 
and 1, all enable tighter determination of the derived kinetic 
constants with varying temperature and concentration. In order 
to achieve such a complex validation technique, additional pro- 
gramming was performed in Lab VIEW. Fig. | shows the fit of 
model to the experimental data from a cell operating at 30, 60 
and 90°C respectively with varying methanol concentrations of 
0.25, 0.5, 1, 2 and 4M. 
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Increasing 
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Fig. 1. Model (line) vs. experimental fit for cell at (A) 30°C, (B) 60°C and (C) 
90°C with varying methanol concentrations. 


The onset potential of the methanol oxidation reaction 
occurred at ca. 0.250 V versus DHE. From Fig. | it can bee seen 
that the kinetic region is not a function of concentration whilst an 
increase in concentration led to increasingly large limiting cur- 
rents that above 2 M were outside the range of the experimental 
data, due to the potential limitation of ruthenium dissolution 
above 0.750 V. The model is based on kinetics and it is essen- 
tial that it fits the experimental data at the low current end of 
the polarisation. Hence, along with the anode polarisation curve 


0.8 


Cell Potential /V 
oO 
A 


0.3 
0.24 
Increasing Temperature 
0.1 
0 
0.1 1 10 100 1000 


Current Density / mAcm-2 


Fig. 2. Model (line) vs. experimental fit on a Tafel plot at 0.25 M (black marker) 
and 0.5 M (white marker) at varying temperature of 30, 60 and 90°C. 


in this work the model was also fitted on a Tafel plot. Fig. 2 
shows the fit of the polarisation data of Fig. 1 for 0.25M and 
0.5 M at varying cell temperature on a Tafel plot which magnifies 
the low current dependence upon potential. The model predicted 
the Tafel slope for the reaction accurately, confirming its reliable 
prediction of the kinetic behaviour. Since experimental factors 
such a membrane electrode manufacturing technique and cell 
test design can affect the mass transport region significantly, in 
this work more attention was given on fitting the data at low end 
of polarisation rather than the mass transport region, which will 
be considered in future publications. 

A further confirmation that the model is functioning correctly 
was the correlation of the kinetic currents of the predicted Tafel 
slope at one potential and three temperatures. The activation 
energy was derived from the Arrhenius relationship and used to 
gauge whether or not the solution was within a suitable range 
between 30 and 70 kJ mol™! [11,22]. Fig. 3 shows the Arrhenius 
plot along with slope for the data generated from Fig. 2. The 
activation energy generated from model was between 42 and 
45 kJ mol~!. Table 1 shows the constants and kinetic parameters 
used in the model. 


5 
y = -5.1703x + 18.858 
4.5] Slope= E, /R 
44 E,=- 42.945 kJ mol 
ee CO Nis 
3} 
<= i 
= 25-4 \ 
Sage 


y = -5.2376x + 18.876 
24 Slope= E,/R 


E= - 43.512 kJ mol’ 


2.5 2.7 2.9 Ai 3.3 3.5 
10° TIK’ 


Fig. 3. Arrhenius plot from the model vs. experimental fit for cell at 0.25 M 
(black marker) and 0.5 M (white marker). 
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Table 1 

Constants and parameters used in the detail kinetic model 

Parameters Fit at 30°C Fit at 60°C Fit at 90°C 

ky 2.6 x 1076 (ms!) 3.3 x 1076 (ms~!) 9.95 x 10-6 (ms!) 

ki 1 x 107!! (mol m~? s) 1 x 107!! (mol m~? s) 1 x 107!! (mol m~? s) 
ky 1.91 x 1078 (mol m~? s) 7 x 1078 (mol m7? s) 1.35 x 1076 (mol m~? s) 
ky 1 x 107! (mol m~? s) 1 x 107! (mol m~? s) 1 x 107! (mol m~? s) 
k3 Ru 6.12 x 1078 (mol m~? s) 5.3 x 1077 (mol m~? s) 4.8 x 1076 (mol m~? s) 
k3 Ru 3.63 x 107! (mol m7? s) 9 x 107! (mol m7? s) 9 x 107! (mol m~? s) 

k4 3.91 x 107? (mol m7? s) 7.62 x 107? (mol m~? s) 9.85 x 107? (mol m~? s) 
a2 0.57 0.65 0.80 

Bo 0.5 0.5 0.5 

B3 0.5 0.5 0.5 

Ba 0.5 0.5 0.5 


4.2. Surface coverage 


The advantage of using such a detail kinetic model lies in the 
determination of surface coverage of intermediate species on 
the catalyst site. Fig. 4 depicts the surface coverage of methanol 
(6m,pt) and hydroxyl ion (@0H.Ru) on Pt—Ru catalyst sites at dif- 
ferent temperature. From these figures, it can be seen that the 
surface coverage of methanol decreases with the increase in cell 


Surface Coverage 


0 "e y 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


Surface Coverage 


0.5 
Potential / V 


° 
N 
a 


Surface Coverage 
o 
oa 


0.4 
Potential / V 


Fig. 4. Surface coverage of intermediate species on Pt—Ru catalyst sites pre- 
dicted by the model with varying methanol concentrations (+: 0.25 M, *: 0.5 M, 
A: 1M, O: 2M, V: 4M) at (A) 30°C, (B) 60°C and (C) 90°C, respectively. 


potential whereas the surface coverage of hydroxyl ion increases 
with the increase in cell potential. It can also be seen that the 
temperature and concentration have opposite effects on hydroxyl 
and methanol coverage on the catalyst site. An increase in tem- 
perature causes 64 pr to decrease and 60H Ru to increase whereas 
an increase in methanol feed concentration causes an increase 
in Oypr and decrease in OoH Ru. Overall, on the basis of the 
response of the surface coverage’s of intermediate species to 
change in temperature and concentration, it can be concluded 
that the anode polarisation curve shows a strong dependence 
on these two operating parameters [5,12,13,17—19,22,23]. From 
Figs. 1 and 4 it can be observed that the limiting current of 
the anode polarisation curve is reached when either the sur- 
face coverage of methanol approaches zero or when the surface 
coverage of hydroxy] group approaches a saturation limit. This 
limiting current can be clearly observed in case of 0.25 and 
0.5 M methanol concentration. The kinetic model also depicted 
the surface coverage of CO on the catalyst site, but over the entire 
potential range Oco pr was close to zero. This is reasonable as the 
1:1 a/o of the Pt:Ru catalyst is commonly employed to ensure 
that there is a high 0oH,Ru presence to react with adsorbed CO at 
source rather than requiring it to diffuse across the crystal sur- 
face and poison the catalyst, as it does on pure platinum. These 
results are in agreement with other literature [11—14,22,24]. 


4.3. Limiting reaction 


It was observed that the variation in the kinetic parameters 
presented in Table | had a significant effect upon the shape and 
fit of the modelled curve. Variation of all of the rate constants 
resulted in an onset potential of ca. 0.250 V versus DHE or above 
in all cases which is in good agreement with the experimental 
data. Fig. 5 shows the effect of varying the most significant 
factors a2, kı and kz on the anode polarisation curve. Varia- 
tion in kı was found only to affect the limiting current of the 
curve as shown in Fig. 5B. Variation in a2 and kz had a signif- 
icant effect upon the onset potential and the Tafel slope of the 
reaction as shown in Fig. 5A and C, respectively. The kinetic 
parameters of the reverse reactions, particularly kį and k4, did 
not affect the polarisation curve or surface coverage of interme- 
diates, although Ky Ru had a minor effect on the polarisation curve 
and surface coverage of intermediates. Rate constants k3 Ru and 
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k4 were found to have a less significant effect upon the curve 
in the fitting range, indicating that, for this model, methanol 
adsorption and oxidation to CO was the rate limiting process 
and that sufficient adsorbed hydroxyl species were available for 
a much faster complete oxidation of CO to CO2. This was also 
demonstrated by low surface coverage of CO on the catalyst site. 


or almost no effect on the polarisation curve. Particularly the 
reverse reaction constants ki and k5, suggest that they can be 
assumed to be negligible in comparison to forward reaction con- 
stants. Exclusion of these constants resulted in a mechanism 
similar to that proposed by Gasteiger et al. [10,14], which sim- 
plified Eq. (6) further to: 


(I-82) FE an FE 


d0 0 0 ¢ ) ¢ ) 
T q Cm (1 Sot 9co,Pt D SO +$efo,ne RT -k Oye PT 


5. Simplified model 


The kinetic model described above fits the experimental data 
well but it is complex, taking a considerable amount of time 
to compute. Variations in the kinetic parameters presented in 
Table 1 depicted that some kinetic parameters have less effect 


Increasing ©, 


Cell Potential / V 


(0) 50 100 150 200 250 300 350 400 
Current Density / mA cm-2 


Increasing k, 


Cell Potential / V 


0 100 200 300 400 500 
Current Density / mA cm-2 


Cell Potential / V 


0 50 100 150 200 250 300 350 400 


Current Density / mA cm-2 


Fig. 5. Effect of varying the influential parameters (A) a2, (B) kı and (C) ko, 
respectively, on the anode polarisation curve. 


At steady state this gives 


ki Cm(1 — 8co,pr) 
ky e((@2FE)/RT) + ky CM 


In the detail kinetic model it was seen that the surface cov- 
erage of CO was very low and is almost negligible above 0.3 V. 
On this basis, it was assumed that, 1 — Oco © 1. Thus, the above 
equation further simplifies to: 


kiCm 
m = ky e((@2FE)/RT) + kı CM 


(5) 
Similarly Eq. (7) further simplifies to: 


picor = kaby e @2FÐ/RT) _ (((1—f4)FE)/RT) 


d k400H, Ruĝ cO, Pt € 


+, kabon ruĝeco pr elC BOFEVRT) — poy el(@2FEV/RT) q6) 


Substituting Eq. (16) into Eq. (14) results in 
z. J = nFkz0M el02FB)/RT) 


Substituting, 0m, from Eq. (15) in above equation and simplify- 


ing it: 
kiCm 
ko e((@2FE)/RT) + k = (17) 


J= nra deronn ( 


Fig. 6 demonstrates the fit of the above simplified model to 
the experimental data of cell at 30, 60 and 90°C with varying 
methanol concentrations. From these figures it can be seen that 
the simplified model also fits the data well but with less accuracy, 
R? statistics 85-90%, as compared to the detailed kinetic model 
with R? statistics 90-95%. 

The model was also fitted to the Tafel plot and the activation 
energy derived in this case were between 42 and 49 kJ mol~!, 
which was slightly higher than the detailed model but still within 
the acceptable range [11,22]. The drawback of the simplified 
model was that it cannot depict the surface coverage of inter- 
mediates involved in the mechanism. The kinetic parameters 
derived from the simplified model are shown in Table 2. 

A further rearrangement of constant kı in Eq. (17) gave: 


n Fky e(@2FE)/RT) CM 
~ Cut (kz / ke (A —2)FE)/RT) 


ae | (18) 
The above Eq. (18) is similar to Meyers and Newman [25] but the 
constants have different interpretations and significance. More- 
over, in this work the above simplified equation is derived step 
by step from the complex kinetic mechanism on the basis of 
proven assumptions. 
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Table 2 
Constants and parameters used in the simplified model 
Parameters Fit at 30°C Fit at 60°C Fit at 90°C 
ky 2.6 x 1077 (ms~!) 3.3 x 10-7 (ms~!) 9.95 x 10-7 (ms7!) 
ko 1.3 x 107? (mol m~? s) 5.04 x 107° (mol m~? s) 1.53 x 1077 (mol m~? s) 
a2 0.60 0.67 0.80 
0.8 ; ; — = o 
0.7 -3 =-1.72x + 2.46 
0.6 4 -6 
2 
S 05 J — -9 
5 2 
o i E 
a 0.4} Increasing 4 42 
3 Concentration y =-2.41x- 5.07 
(S) WORMS See eS es seq uses 
os TTT Fha ee eee ee 15) D SNN seen O y =-2.41x- 7.38 
MES Se, y =-7.97x + 9.64 
0.2 “8p 
A oY =-8.60x + 7.56 
Oi i ne eR -21 - - <= 
0 50 100 150 200 250 Zr 2.9 3.1 3.3 3:5: 
Current Density / mAcm-2 10TIK 


Increasing 
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Cell Potential /V 


0 100 200 300 400 500 600 
Current Density / mAcm-2 


Increasing 
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Current Density / mAcm-2 


Fig. 6. Model (line) vs. experimental fit for cell at (A) 30°C, (B) 60°C and (C) 
90°C with varying methanol concentrations. 


Fig. 7. Arrhenius plot of the kinetic parameters for comparison of resulting 
activation energy with other published results (@: ko, W: k3 Ru, A: k4, O: ky 
(simplified model), O: k2 (simplified model)). 


Finally the parameters obtained by both the detailed model 
and simplified model were compared with other published 
results. The resulting kinetic parameters from both the modelling 
procedure were in accord with other literature [12,13,25,26]. 
Moreover, the activation energy obtained by plotting the Arrhe- 
nius plots, as shown in Fig. 7, also agreed with the published 
results [11,19,22]. Based on the previous work [17,18] and the 
models derived here it was found that the simplified model, 
given by Eqs. (17 and 18), involved the three significant kinetic 
parameters (a2, kų and k2) which were mainly responsible 
for describing the behaviour of anode polarisation curve as a 
function of cell temperature, anode potential and methanol con- 
centration. The detailed models, given by Eq. (14) and Eq. (A.2), 
due to the presence of above significant parameters describe the 
polarisation behaviour and in addition, due to the presence of 
other kinetic parameters, it’s useful in revealing the surface cov- 
erage intermediates on the catalyst site. The detailed models 
can give slightly more accurate results over the full range of 
anode potentials, although the simpler model is sufficient for 
most simulations. 


6. Conclusion 


A complex methanol oxidation mechanism was studied with 
the help of detail kinetic model. The surface coverage of inter- 
mediate species formed during the reaction were calculated and 
agreed with the results from other literature [12,13,19]. The 
model calculates anode polarisation behaviour on the basis of 
surface coverage and shows that the limiting current is reached 
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when the surface coverage of adsorbed methanol approaches 
zero or when the surface coverage of hydroxy] group approaches 
the saturation limit. The model highlights the subtle balance 
between the methanol adsorption and subsequent oxidisation 
to CO on the dual site catalyst and identifies them as the lim- 
iting steps rather than the surface oxidation of CO to CO2. It 
was observed that by varying the kinetic rate constants a surface 
coverage of intermediates species can be manipulated while still 


oj = nF (ka POOH, pi + k4, Ru90H,Ru)Oco et PH FER 


(A.2) 


By assuming the activity of water equal to unity the steady state 
solution in this case can be further simplified to: 


ky Cu — Om(kiC + kj + kz ef(@2FE)/RP)) 


—k,Cm@ou,Pt 


obtaining a reasonable fit to the anode polarisation curve. To fur- OCO, Pt = LG Ee ADFERT) (A.3) 
ther validate the kinetic constants and the modelled coverage at IM 82% 
particular operating conditions, the behaviour of intermediate ko ef (02FE)/RT) 
species with respect to temperature, concentration and potential Gco,P = k, e(1—B2)FE)/RT) 4 ky Oo pt e(U—B4)FE)/RT) = 
should be confirmed by the use of modern spectroscopy tools like + kabon gu eC! BH)FED IRT) 
7 k3.p(1 — Oon.pt — Opel B3)FE)/ RT) = kh pOoH.Pt e(((—B3)FE)/ RT) " 
CO,Pt = k3 pi e(-B3)FE)/RT) + k4 pOoyy pt e(((l-BaFE)/RT) a 
k3 Ru(l — OOH, Ru Jel Bs) FE)/ RT) = k% RuÎOH, Ru e(((—B3)FE)/RT) 
9co,Pt = (A.6) 


the Fourier Transform Infra-red Spectroscopy (FTIRS) [3,20,24] 
and corroborated with the model. 

Further refinement of this model can be achieved by adding 
more realistic conditions to the original reaction expressions. 
This could include non-ideality of the platinum and ruthenium 
distribution which would influence the impact of k3 and k4 
and the addition of rate expressions for the more significant 
complex and stable intermediates. This aspect of modelling is 
the subject of ongoing work in the laboratories at Newcastle 
University. 

The simplified model derived from the complex mechanism 
gives a faster approximation of anode performance. This would 
be useful if coupled with a diffusion model for the prediction of 
performance along a flow field channel or through the Nafion® 
membrane. Overall, depending on the computational speed and 
complexity of problem at hand, either of the models can be used 
for complete cell or stack studies using kinetic data. 
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Appendix A 


In this modelling procedure the surface coverage of hydroxyl 
ion on Pt site is considered and it is assumed that the reaction of 
OHads With COaas,pt on both Pt and Ru sites occurs at different 
rates which adds another two kinetic parameter to the model. 
Thus, the rate of reaction is given by: 
r4 = ka,pOou,prOco el PUFED/RT) 


+ ka RuP0H Rub co eM Pa FEV/RD) (A.1) 


kalok Ru eI —BaFE)/RT) 


These non-linear equations can be solved further by using 
mathematical software like Maple and the solution can be 
incorporated into (A.1) and (A.2) to find the respective rate 
constant and subsequently the current density [17]. This mod- 
elling procedure will contain additional three kinetic parameters 
(k3,Pt, k$ pp k4pt) hence it will be more time consuming and 
complex compared to the detail model in manuscript. 
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